The wave function of [70, 1 ] baryons in the l/N^ expansion 



N. Matagne* and Fl. Stancu^ 
University of Liege, Institute of Physics B5, 
Sart Tilman, B-^OOO Liege 1, Belgium 
(Dated: February 26, 2008) 

Abstract 

Much work has been devoted to the study of nonstrange baryons belonging to the [70, 1~] 
multiplet in the framework of the 1/Nc expansion. Using group theoretical arguments here we 
examine the relation between the exact wave function and the approximate one, customarily used 
in applications where the system is separated into a ground state core and an excited quark. We 
show that the exact and approximate wave functions globally give similar results for all of mass 
operators presented in this work. However we find that the inclusion of operators acting separately 
on the core and on the excited quark deteriorates the fit and leads to unsatisfactory values for the 
coefficients which encode the quark dynamics. Much better results are obtained when we include 
operators acting on the whole system, both for the exact and the approximate wave function. 
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I. INTRODUCTION 



The 1/Nc expansion of QCD 



approach to study 

baryon spectroscopy. It has been apphed to the ground state baryons QM y, 0, y, y, q 

as well as to excited, states, in particular to the negative parity [70, 1~] multiplet (A^ = 1 
band) [n, 12, 3, 



14 llSl . llGl]. to the positive parity Roper resonanc e [5 6', 0^] (A^ = 2 band) 
17|, to the [56,2+] (iSj and the [70, £+] multiplets (£ = and 2) [igi |2o|, both belonging 



to the N = 2 band and to the [56,4+] multiplet (A^ = 4 band ) |2l[| . Estimates for the 
lowest multiplet [70, 3^] of the = 3 band have also been made 22 [.In this approach the 
main features of the constituent quark model emerge naturally |23|, l2J] as for example the 
dominant role of the spin-spin term and the smallness of the spin-orbit term. 

The study of excited states belonging to the symmetric representation [56] is similar to 
that of the ground state. The introduction of an orbital part in the wave function does 
not affect the procedure, inasmuch as the flavor-spin part remains symmetric. The study 
of excited states which are mixed symmetric, both in the orbital and flavor-spin space is 
more complicated and it became controversial. Presently we discuss the [70, 1~] multiplet, 
which is the simplest one, having 7 experimentally known resonances 25|] in the nonstrange 
sector. The standard procedure 13j is to separate the system into of a ground state core 
of A"c — 1 quarks and an orbitally excited {£ = 1) quark. Then, in the spirit of a Hartree 
picture the system is described by an approximate wave function, where the orbital part 
has a configuration of type s'^'^'^p (no antisymmetrization) which is combined with an 
approximate spin- flavor part. There is also a more straightforward procedure where the 



system is treated whole 



261 ]. This procedure is in the spirit of the spectrum generating 



algebra method, introduced by Gell-Mann and Ne'eman. It has been applied to the ground 
state, where use of the generators describing the entire system is made. In this procedure 
the exact wave function is needed. 

Here we reanalyze the [70, 1~] multiplet (A^ = 1 band) without any prejudice. We start 
with the standard procedure based on the core+excited quark separation. We show the 
relation between the approximate wave function and the exact one, constructed explicitly in 
the next section. To test the validity of the approximate wave function 13| we first compare 
the analytical expressions of the matrix elements of various operators entering the mass 
formula. Next we perform a numerical fit to the data. For practical reasons we choose a 
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simple mass operator containing the most dominant terms in order to fit the 7 experimentally 
known masses and ultimately the two mixing angles also. 

The next section is devoted to the analysis of the wave function and its construction by 
using isoscalar factors. In Sec. Ill we briefly introduce the mass operator. In Sec. IV we 
perform fits to the data, by including various sets of invariant operators constructed from 
the generators of SU(4) and 0(3). Sec. V is devoted to mixing angles. The last section 
contains our conclusions. Appendix A describes the procedure used to construct the needed 
isoscalar factors of the permutation group. Appendix B introduces the fractional parentage 
coefficients needed to separate the orbital wave function into a part describing the core 
formed of — I quarks and a single quark. Appendix C provides analytic expressions of 
the matrix elements necessary to calculate the isoscalar factors. 

II. THE EXACT WAVE FUNCTION 

The system under concern has four degrees of freedom: orbital (O), flavor (F), spin (S) 
and color (C). It is then useful to construct its wave function with the help of inner products 
of in order to fulfill the Fermi statistics. The color part being antisymmetric, the orbital- 
spin-fiavor part must be symmetric. As in the mass formula there are no color operators, 
the color part being integrated out, we are concerned with the orbital-spin-fiavor part only. 

A convenient way to construct a wave function of given symmetry [/] is to use the inner 
product [/'] X [/"] of two irreducible representations (irreps) of S^^, generally reducible into 
a Clebsch-Gordan series with a given number of irreps [/]. A basis vector of [/] is denoted 
by I [fjy), where Y is the corresponding Young tableau (or Yamanouchi symbol). Its general 



where S{[f']Y'[f"]Y"\[f]Y) are Clebsch-Gordan (CG) coefficients of Sn,. 

Here we deal with a system of Nc quarks having one unit of orbital excitation. Hence 
the orbital part must have a mixed symmetry [Nc — 1, 1]. To get a totally symmetric state 
[[Ai'c]!), where for brevity the Yamanouchi symbol {Nc times) is denoted by 1, the FS 
part must belong to the irrep [Nc — 1, 1] as well. Then the wave function takes the form 



form is 27 1 



\[f]Y)= Y: S{[f']Yy'']Y''\[f]Y)\[f']Y%ny''), 



(1) 



Y',Y" 




1/2 



(2) 



Y 
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where the coefficient in front is the CG coefficient needed to construct a symmetric state of 
Nc quarks. The sum is performed over all possible Nc — 1 standard Young tableaux. In this 
sum there is only one Y (the normal Young tableau) where the last particle is in the second 
row and Nc — 2 terms with the A^^c-th particle in the ffist row. If we specify the row p of the 
Nc-th quark and the row q of the {N^ — l)-th quark, and denote by y the distribution of the 
Nc — 'i remaining quarks, we can write Y more explicitly as 

Y = (pqy). (3) 

For illustration let us consider the case A^^^ = 5. Expressed in terms of Young tableaux, the 
wave function (El) reads 
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(4) 



In this notation, in the left hand side one has p = q = 1 because the A^^c-th particle, i.e. 
particle 5, and the {Nc — l)-th, i.e. particle 4, are both in the ffist row. The ffist term in 
the right hand side of (HI) has p = 2 and q = 1 both for (O) and (FS) parts. This is the only 



term considered in Ref. 13|. The second term has p = 1 and q = 2 and the last two terms 
have p = q = 1. In the following we shall mostly concentrate on the p = 1 terms in the right 
hand side, neglected in Ref. [l3| and compare results derived from the exact wave function 
([2]) with those obtained in Ref. [l^. To include the neglected terms we have to explicitly 
consider the (FS) part of the wave function, as follows. 

At its turn, the (FS) part can be decomposed into a sum of products of spin (S) and flavor 
(F) parts using Eq. ([1]) where in the left hand side we now have = |[A'c — 1, 1]Y)fs 

and in the right hand side the appropriate CG coefficients. For these coefficients it is useful 
to use a factorization property (Racah's factorization lemma) which helps to decouple the 
A"c-th particle from the rest. Accordingly, every CG coefficient of Sn^ can be factorized into 
an isoscalar factor times a CG coefficient of Sn^-i, and so on. We apply this property to 
3 'S'at^^i only. To do so, it is necessary to specify the row p of the A'c-th and the row q 
of the (A'c — l)-th quarks, as above. If the isoscalar factor is denoted by K{[f']p'[f"]p"\[f]p) 
the factorization property reads 



where the second factor in the right hand side is a CG coefficient of Sat^-i containing the 
partitions [ fL] , [fp„] and [fp] obtained after the removal of the Ai'c-th quark. 



In Ref. [19(] we have derived the isoscalar factors associated to the normal Young tableau 
of \[Nc — 1,1]Y)psj "i-G- for P = 2, g = 1. We have also shown that the coefficients Cp^ 
appearing in the wave function of Ref. [3] correspond precisely to these isoscalar factors. 
In the present work we derive the isoscalar factors of the other Nc — 2 terms, having p = 1, 
for the ffist time. Details are given in Appendices A and C. The results are summarized 
in Tables I, II and III. In each table the column corresponding to p = 1 is new. For 
completeness, we also reproduce the results of Ref. isl for p = 2. They are both used 
in the calculation of matrix elements of various operators entering the mass formula. As 
shown below, the total wave function needed to calculate these matrix elements contains a 
spin-flavor part, Eq. ([8]), constructed with the help of the isoscalar factors of Tables I, II 
and III. 

The isoscalar factors derived here can be applied to the study of other [70, i^] multiplets 
with i 1 and parity P = + or P = — , or to other physical systems having fermions. 



They satisfy a recurrence relation described in Refs. |27l . l28l | which has been used to check 
the presently derived analytic results in the case of A",, = 5, 7 and 9 quarks. The results 
from Tables I, II and III can also be interpreted as isoscalar factors of SU(4), related to CG 
coefficients appearing in the chain SU(4) D SUs(2) x SU/(2) with the specific values of S 
and / presented here. 

For the orbital part of the wave function ([21) one can use one-body fractional parentage 
coefficients (cfp) to decouple the A^c-th quark from the rest. Together with the factorization 
property ([5]) it leads to a form of the wave function which is a sum of products of a wave 
function describing a system of A^c — 1 quarks, called "core" times the wave function of the 
A'c-th quark. If couplings are introduced this gives 

\iSJJs;Ih) = E ci{p,ic,Q \WSJJ3;Ih)p, (6) 

where a{p, £c, iq) are the one-body cfp given in Appendix B, with 4 and ig representing the 
angular momentum of the core and of the decoupled quark respectively. One can write 



\WSJJs-Ih)p 
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\icm,)\igmg)\[N, - 1, l]p- SSs; II3), (7) 



containing the spin-flavor part 



|[iV, - 1, l]p; SSs; Ih) = E KmPify'liNc - 1, l]p)\SSs;p')\Ih;p") 



(8) 



where 



\SSs;p')= E 



nil 



Sr. 



S 



\Scmi)\l/2m2), 



(9) 



with Sc = S- 1/2 if p' = 1 and Sc = S + 1/2 if p = 2 and 

|//3;/) = E 2 ' 1 |/,Z,)|1/2Z2), 

n,i2 \ ii 12 h 



(10) 



with Ic = I — 1/2 if p" = 1 and Ic = 1 + 1/2 if p" = 2. Thus p' and p" represent 
here the position of the Nc-th quark in the spin and isospin parts of the wave function 
15* 5*3; p')\IIs; p") respectively. Note that the wave function (jS]) should also contain the second 
factor S{\f'pi]q' y'[fpi]q" y"\[fp]qy) of Eq. ^ together with sums over q'y' and q"y" which 
ensures the proper permutation symmetry of the spin-flavor part of the core wave function 
after the removal of the A^^c-th quark. However, these CG coefficients of Sat^-i do not need 
to be known and for simplicity we ignore them because their squares add up together to 
1, by the standard normalization of CG coefficients. Equations fl6|)- f|T0l) are needed in the 
calculation of operators as s ■ Sc, S^, t ■ Tc, T^, etc. (see next section) which involve the core 
explicitly. 

Let us illustrate, for example, the use of Table I, when p = 1. Then Eq. ([8]) becomes 
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(11) 



where Sc and Ic are determined as explained below Eqs. © and (ITU]) . 

In the following we shall compare results obtained from the exact wave function ^ with 
those obtained from a truncated wave function, representing only the term with p = 2 in 



the sum ([6]). In addition, by taking a(2,£c = 



1 



1, instead of the value given in 



Appendix B, and renormalizing the resulting wave function, we recover the wave function of 
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TABLE I: Isoscalar factors K{[f']p'[f"]p"\[f]p) for S = I = 1/2, corresponding to ^8 when Nc = 3. 
The second column gives results for p = 1 and the third for p = 2. 
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TABLE IL Isoscalar factors K(\f'\p'[f"\p"\[f\p) for S = 3/2, / = 1/2, corresponding to ^8 when 
Nc = 3. The second column gives results for p = 1 and the third for p = 2. 



Ref. ^3|]. As seen from Appendix B, it is entirely reasonable to take a(2, £c = 



when Nc is large. Then the core is in its ground state. It is in fact the validity of this wave 
function, representing the Hartree approximation, we wish to analyze as compared to the 
exact wave function in the description of the [70, 1~] multiplet. 



7 



[fWV 



[iVc-i,i]i 



[A^c-l,l]2 





Nc- 


1' 


1 


riVc + 3 


iVc- 


3' 


1 




/(iV,-l)(iVe + 3) 


[ 2 ' 


2 




[ 2 ' 


2 






/ Nc{Nc - 2) 




Nc- 


!■ 


2 


riVc + 3 




3' 


2 




/5(iVe-l)(iVe-3) 


[ 2 ' 


2 




[ 2 ' 


2 






2N,{N, - 2) 






1' 


2 


\N, + 3 




3' 


1 




/(iV,-3)(iV, + 3) 


[ 2 ' 


2 




[ 2 ' 


2 






/ 2Nc{Nc - 2) 



2 ' 2 

TABLE III: Isoscalar factors K{[f']p'[f"]p"\[f]p) for S = 1/2, / = 3/2, corresponding to ^10 when 
Nc = 3. The second column gives results for p = 1 and the third for p = 2. 

III. THE MASS OPERATOR 

The mass operator M is a linear combination of independent operators Oj 

M = 5^cA, (12) 

i 

where the coefficients Cj are reduced matrix elements that encode the QCD dynamics and 
are determined from a fit to the existing data. The operators Oi are constructed from the 
SU(4) generators Si, Ta and Gia and the S0(3) generators fj. For the purpose of the present 
analysis it is enough to restrict the choice of the operators Oi to a few selected dominant 
operators. Our previous analysis j:2fi] suggests that the dominant operators up to order 
0{1/N^ included are those constructed from the SU(4) generators exclusively, the operators 
containing angular momentum having a minor role. Among them, the spin-orbit, although 
weak, is however considered here. Samples are given in Tables VI-X. (For convenience and 
clarity we denote some operators in the sum (fT2l) by 0[. The meaning will be obvious later 
on.) 

In Tables IV and V we present the analytic expressions of the matrix elements of the 
operators containing spin and isospin, used in this study. They have been obtained either 
with the approximate or with the exact wave function. One can see that the analytic forms 
are different in the two cases. This implies that at A^c = 3 the values obtained for q from 
the fit to data are expected to be also different, as illustrated in the next section. The 
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TABLE IV: Matrix elements of the spin operators calculated with the approximate and the exact 
wave functions. 
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TABLE V: Matrix elements of the isospin operators calculated with the approximate and the exact 
wave functions. 



matrix elements of the spin-orbit operator are i 
wave function and can be found in Table II of Ref. 



dentical for the exact and an approximate 

3- 



IV. FIT AND DISCUSSION 



In the numerical fit we have considered the masses of the seven experimentally known 
nonstrange resonances belonging to the [70, 1~] multiplet. Neglecting the mixing 
we identify them as: 2iVi/2(1538 ± 18), ^A^i/2(1660 ± 20), ^N:i/2{lh2?> ± 8), ^Ar3/2(1700 ± 50), 
^Ar5/2(1678 ± 8), 2Ai/2(1645 ± 30) and 2A3/2(1720 ± 50). 

The results from various fits are presented in Tables VI-IX. In Table VI the first four 
operators are among those considered in Ref. [isl, where the system is decoupled into a 
ground state core and an excited quark. To them we have added the isospin operator which 
we expect to play an important role. The value of Xdof satisfactory but the value of the 
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coefficient ci is mucli smaller than in previous studies where it generally appears to be of 

the order of 500 MeV. The spin-orbit coefficient is small, as expected, but the coefficients 

C3 and C4 are exceedingly large in absolute values and have opposite signs, which suggest 

some compensation. The coefficient C5 is also very large and negative for the approximate 

wave function but it has a reasonable value for the exact wave function. A solution is 

either to eliminate the isospin operator, but this cannot be theoretically justified or to make 

some combinations of the above operators. In Table VII we consider the linear combination 
• 3 

2s^Sl + SlSl + - = 5"^, i.e. we use the total spin operator of the system, entirely legitimate 

in the 1/Nc expansion approach. The coefficient ci acquires the expected value and the 

spin coefficient Cg has a reasonable value, being about 70 MeV smaller for the exact wave 

function than for the approximate one. The isospin coefficient is equal in both cases and 

also has a reasonable value. By analogy, in Table VIII we include the linear combination 
3 

2t"'T^ _l_ T^T^ _|_ _ — /2 ^ g ^gg ^jig total isospin and restrict the spin contribution to 

—s ■ Sc- Then the isospin coefficient Cg acquires values comparable to the spin contribution 
in Table Vll. The spin coefficient C3 has now values comparable to the isospin coefficient 
C5 in Table VII. One can infer that the spin and isospin play a similar role in the mass 
formula. This statement is clearly proved in Table IX where we include the spin and isospin 
contributions on an equal footing, i.e. we include operators proportional to S'^ and P. The 
situation is entirely identical for the exact and approximate wave function, with reasonable 
values for all coefficients and a Xdof = 1-04:. The identity of the results is natural because 
both the approximate and the exact wave function describe a system of a given spin S and 
isospin /. 

V. MIXING ANGLES 

So far we have discussed the mass spectrum. Additional empirical information come from 
the mixing angles extracted from the electromagnetic and strong decays, in particular from 
the dominance of Nr] decay of A^i/2(1538 ± 18). These angles are defined as 

lA^j(upper)) = cosOjI^Nj) + sinOjl'^Nj), 

I iVj (lower)) = cos OjI'^Nj) - sin 6 j\'^Nj). (13) 
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Operator 




Approx. w.f. (MeV) 


Exact w.f. (MeV) 


Ol=iVel 


Ci = 


211 ±23 


299 ± 20 


O2 = ts' 


C2 = 


3± 15 


3 ±15 


1 

O4 = — S^S^ 


C3 = 
C4 = 

C5 = 


-1486 ± 141 

1182 ± 74 
-1508 ± 149 


-1096 ± 125 
1545 ± 122 
417 ± 79 


Xdof 




1.56 


1.56 



TABLE VI: List of operators and coefficients obtained in the numerical fit to the 7 known exper- 
imental masses of the lowest negative parity resonances (see text). For the operators O3, O4 and 
O5 we use the matrix elements from Tables IV and V. 



Operator 




Approx. w.f. (MeV) 


Exact w.f. (MeV) 


Oi = N,l 


Cl = 


513 ±4 


519 ±5 


O2 = ts' 


C2 = 


3 ±15 


3 ±15 




4 = 

C5 = 


219 ± 19 
417 ± 80 


150 ± 11 
417 ±80 


Xdof 




1.04 


1.04 



TABLE VII: Same as Table VI but for O3, which combines O3 and O4 instead of using them 
separately. 



Operator 




Approx. w.f. (MeV) 


Exact w.f. (MeV) 




Cl = 


516 ±3 


522 ±3 


O2 = ts' 


C2 = 


3± 15 


3 ±15 




C3 = 

4 = 


450 ± 33 
214 ± 28 


450 ± 33 
139 ± 27 


Xdof 




1.04 


1.04 



TABLE VIII: Same as Table VII but combining isospin operators instead of spin operators. 
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Operator Approx. w.f. (MeV) Exact w.f. (MeV) 

01 = Ncl ci = 484 ±4 484 ±4 

02 = i's' C2= 3 ± 15 3 ± 15 
0', = ^^(^2s'Si + SlSl + ^^ 4= 150 ±11 150 ±11 

oj, = ^ ^2t"r^" ± r^"r^" ± ^ 4 = 139 ± 27 139 ± 27 

Xdof 1-04 1.04 

TABLE IX: Same as Table VI but with operators proportional to the SU(2)-spin and SU(2)-isospin 
Casimir operators (see text). 

rn 

Experimentally one finds ^^Jg ~ —0.56 rad and 6'^^^2 ~ O-^O [29]. In our case, the only 
operator with non-vanishing matrix elements ^A^ — ^A^, thus contributing to the mixing, is 
the spin-orbit operator. Its off-diagonal matrix elements are (^Ni/2\i ■ s\'^Ni/2) = —1/3 and 



( ^3/2!^ ■ s\ N3/2) = — yS/lS, compatible with Ref. [l3j when Nc = 3. In this notation the 
physical masses become 

Mj(upper) = M^Nj) + C2{^Nj\i ■ s\''Nj) tanOj, 

Mj(lower) = M^Nj) - C2{'^Nj\i ■ s\'^Nj) tanOj. (14) 

Due to the fact that the diagonal contribution of O2 = ^ ■ s is very small, we expect its 
contribution to the mixing to be small as well. By including the experimental values of 6j in 
Eqs. ( 11^ and applying again the minimization procedure, we found indeed that the effect 
of mixing is negligible. The only coefficient being slightly affected is C2 and the Xdof remains 
practically the same. Next we varied 9j to see if the value of Xdof remains stable. Indeed it 
does, for variations of 6j of ±0.05 rad. In Table X we show the optimal set of coefficient s a 



associated to the experimental values of dj. In a more extended analysis as that of Ref. 26] 
the operators O5 = ^i(.^)iiG"'G^'' and Og = J-f T^G*" would also contribute. However 
all these operators contain angular momentum while in constituent quark models where the 
predictions are generally good, the source of mixing are either the spin-spin or the tensor 
interactions (no angular momentum), contrary to the present case. We therefore believe 
that this important issue deserves further work in the future, combined with a strong decay 
analysis, like in Ref. [sol . 



12 



Operator Approx. w.f. (MeV) Exact w.f. (MeV) 

01 = Ncl Cl = 484 ±4 484 ±4 

02 = i's' C2= -9 ± 15 -9 ± 15 
0', = ^^(^2s'Si + SlSl + ^^ 4= 150 ±11 150 ±11 

oj, = ^ ^2t"r^" ± r^"r^" ± ^ 4 = 139 ± 27 139 ± 27 

Xdof 1-05 1.05 

TABLE X: Fit with state mixing for and due to the spin-orbit couphng i ■ s. 

VI. CONCLUSIONS 

This study shows that the exact and the approximate wave function give identical Xdof 
in this simphfied fit and rather similar results for the dynamical coefficients q entering in 
the mass formula. It also shows that the contribution of the isospin operator O5 = j^^^ 

is as important as that of the spin operator O3 = j;^^"^- addition, it turned out that 
the separation of S'^ and of into independent parts containing core and excited quark 
operators is undesirable because it seriously deteriorates the fit. The only satisfactory way 
to describe the spectrum is to include in Eq. ( fT2l) the Casimir operators of SU5(2) and 
811/(2) acting on the entire system. 

The difficulty with the approximate wave function is that it cannot distinguish between 
S'^ and J^, which act on the core only. In practice it means that the contribution of j^^c 

the mass is hidden in the coefficient q associated to — 5^. This also happens for — for 
states described by a symmetric flavor-spin wave function where its contribution is absorbed 
by the coefficient of — S'^. In oti 
mass by about the same amount. 

In the decoupling scheme the is 
which manifestly deteriorates the fit, as seen in Table VI. This may explain why this operator 



by the coefficient of — S"^. In other words S*^ and should share their contribution to the 

Nr 

e amount. 

In the decoupling scheme the isospin can be introduced only through the operator —t-Tc 



has been avoided in numerical fits of previous studies based on core ± quark separation [13 1 



While preserving the S^^ symmetry in these calculations, our conclusion is at variance 
with that of Pirjol and Schat [31] who formally claim that the inclusion of core and excited 
quark operators is necessary, as a consequence of constraints on the mass operator, resulting 



13 



from the same S^^ symmetry. 

The dynamical coefficients q encode the quark dynamics in hadrons. Therefore, to find 
their correct values and their evolution with the excitation energy is a very important task. 
They can serve to test the validity of quark models and quantitatively determine the contri- 
bution of gluon and Goldstone boson exchange interactions, both justified by the two-scale 



picture of Manohar and Georgi 32|. 



As a byproduct we have obtained isoscalar factors of the permutation group which may 
be useful in further studies of mixed symmetric multiplets within the 1/Nc expansion or in 
any other approach where one fermion is separated from the rest. 

Acknowledgments. The work of one of us (N. M.) was supported by the Institut 
Interuniver sit aire des Sciences Nucleaires (Belgium). 

APPENDIX A: ISOSCALAR FACTORS 

Here we shortly describe the derivation of the isoscalar factors of Tables I, II and III. We 
ffist recall that the expressions of K{[f']p'[f"]p"\[Nc — 1, l]p) with p = 2 were obtained in 
Ref. jl9|. We are presently concerned with the calculation of K{[f']p'[f"]p"\[Nc — 1, l]p) for 
p = 1. 

In the derivation below we use the symmetry property 

KiifWYimp) = Kiirvimmp), m 

and the orthogonality relation 

E^([/>'[/"KI[/]p)^([/>'[/"]/l[/iK) = ^fh^ppi- (A2) 

p'p" 

1. The case S = I = 1/2 

In this case we have 



[/I = [/" 



(A3) 



2 2 

and there are three nonzero coefficients necessary to express the FS state \[Nc — 1, as 
linear combinations of S and F states. These are K{[f]2[f"]2\[N^-l, 1]1), K{[f]l[f"]2\[N^- 
1, 1]1) and K{[f']2[f"]l\[Nc - 1, 1]1). But due to (jSH) the latter two are equal. Then we are 
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left with two unknown isoscalar factors which following (lA2p are normalized as 

{KiimmN, - 1, + 2iKmnnmc - 1, mr = i. (A4) 

To find the solution we need one more equation. This is provided by the matrix elements 
of the operator G"" written as 

{G'-) = {G';) + (g^-) (A5) 

with G*" and g^"" acting on the core and on the excited quark respectively. According to 
Eqs. f ICip and (IC3I) the matrix elements of the SU(4) irrep [Nc — 1, 1] with S = I = 1/2 are 

G'" -S3; -/3\ = 
2 2 / 

/ 



2 •^3.2^3 



12 



V 



1/2 1/2 1 



1/2 



4 



(A6) 



1/2 1 

Taking p = 1 in Eqs. fICil) and fICSl) one obtains for (G*") and (^f*") the following 
expressions 



Gl 



{Kiif'wnN, - 1, i]i))' + (N, + 3) {Kiimnmc - 1, i]i))' 



+ 4 
/ 

X 

V 



'(AT, -3) (AT, -1) 



Kiimnmc - 1, i]i)K([/']i[r]2|[Ar, - 1, i]i) 



1/2 1 

S3 i 



1/2 W 1/2 1 



S' 



h a 



1/2 

^3 



(A7) 



and 



[N, - 1, 1]1; \s',- h', |r I [Nc - 1, 1]1; \S3; \h 



^ (;^([/']2[r]2|[Ar, - 1, 1]1))^ - i {K{[mf"]2m - 1, 1]1))' 



X 



1/2 1 

5*3 i 



1/2 W 1/2 1 



S' 



L a 



1/2 

^3 



(A8) 



respectively. Then Eq. f lASp leads to the following relation 



+ ^ir([f ]2[r]2|[Ar, - 1, 1]1))2 + ^(K([f ]l[r]2|[Ar, - 1, 1]!)^ 



12 



6 



2 / (Ar,-3)(Ar,-l) 
'3V 2 

A-c-e 

— = 

12 



KiimmN, - 1, i]i)ir([/']i[r]2|[Ar, - 1, i]i) 



(A9) 
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which together with the normahzation relation (lA4p forms a system of two nonhnear equa- 
tions for the unknown isoscalar factors. The solution is exhibited in Table I column 2. The 
phase convention for the isoscalar factors is the same as in Refs. |27l. |28|. 



2. The case S = 3/2, I = 1/2 



These values of S and / imply 



[/'] 



[/"] 



iVc + 1 iV,-l 
2 ' 2 



(AlO) 



and, as above, there are three nonzero coefficients necessary to express the FS state | [N^ — 
1, as linear a combination of S and F states. In this case, these are i^'([/']2[/"]2| [Ai'c — 
1, 1]1), K{[f]l[f"]2\[N^ - 1, 1]1) and K{[f]l[f"]l\[N^ - 1, One needs three equations 
to find the solution. 

First, Eq. (lA2p gives the normalization relation 



(K([r]2[n2|[Ar,-l,l]l))2 + (K([r]l[/n2|[Ar,-l,l]l)f + (K([/']l[r]l|[iV,-l,^ 

(All) 



Second, Eq. (lASP gives 
3{3N^ - 1 



iVc + 1 



{Kiimnm, - 1, i]i))^ + ^—(Kumnm, - 1, i]i)) 



ION, 



2Nr. 





(iVc- 


-l)(^^c + 3)^ 










-3)(iVe-l) 














5N^ 



i^([/']l[r]2|[Ar, - 1, l]l)ir([/']l[r]l|[Ar, - 1, 1]1) 
ir([f]2[r]2|[iV,-l,l]l)ir([/']l[r]l|[iV,-l,l]l) = l.(A12) 
Third, for the irrep [/] = [A^^c — 1, 1] the Casimir operator identity becomes (see Eq. ( ]C3[) ) 



+ 2 (a 



2 2 

from which one can derive the following equation 

{s-S,) + {t-T,)+4{g-G,) 
16 



N,{3N, + 4) 



(A13) 



3N,-7 
4 ■ 



(A14) 



Using Eqs dCSD-dCHD we obtain 



3K 



{Ki[f']2[nmc-i,mf 



iVc + 7 
8 



{K{[f']i[nmc-i,mr 




Kiimnmc - 1, mKumnmc - 1, i]i) 




K{[f']l[f"]2\[N, - 1, l]l)K{[f]l[ni\[Nc - 1, 1]1) 



Kiimnmc - 1, i]i)Ai[/']i[r]i|[iv, - 1, i]i) 



(A15) 



which is the third equation needed to derive the isoscalar factors presented in the second 
column of Table [Tll 

3. The case S = 1/2, I = 3/2 

This case is similar to the previous one. The results can be obtained by interchanging S 
with /. 

APPENDIX B: FRACTIONAL PARENTAGE COEFFICIENTS 

We are in the case where the orbital wave function \[Nc — 1, l]pqy)o contains only one 
excited quark with the structure s^^^^p. The one-body fractional parentage coefficients (cfp) 
help to decouple the system into a core of Nc — 1 quarks and a single quark. 

For p = 2 { obviously g = 1 and y is fixed) the decoupling into a core and a quark gives 



\[Nc - 1, l]2)o = - (Bl) 



In the first term the core is in the ground state, i.e. ic = 0, and is described by the 



momentum iq = 1. In the second term, the orbital wave function is still symmetric but it 
contains a unit of angular momentum and the quark is unexcited. 
For p = 1, irrespective of q and y one has 





symmetric orbital wave function R^j^^_^^{s 



Nc-l 



) . The quark is excited and carries the angular 



\[Nc-l, 



l]l)o - -R[7Ve-2,l]('5 



(B2) 
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i.e. the core is always excited, thus it is described by a mixed symmetric state, and the 
quark is in the ground state. Denoting the one-body cfp by a{p,ic,^q) we get 



to be used in Eq. 



APPENDIX C 



N -1 

a(2,4 = 0,£, = l)= ' ' 



a(2,4 = 1,^, = 0) = - 



a{l,i, = ljg = 0) = 1. 




(B3) 

(B4) 
(B5) 



This appendix contains the matrix elements of some operators used to derive the isoscalar 
factors K{[f']p'[f"]p"\[N^ - 1, 1]1) presented in Tables [H [TTl and [ml 

Let us first recall the general formula of the matrix elements of the generator G*" of SU(4) 
for a given irrep [/]. This is [33 1 



{[f];I'l',;S'S',\G'''\[f];Ih;SSs) 
/CI/](SU(4)) ( [f] [21' 
IS 11 




p=i 



S 1 

S3 Sn 



S' 
S' 



I 1 



I' 



(CI) 



where C^'^' is the eigenvalue of the SU(4) Casimir operator for the irrep [/], the second factor 
is an isoscalar factor of SU(4), the third a CG coefficient of SU(2)-spin and the last a CG 
coefficient of SU(2)-isospin. The isoscalar factors used here are those of Table A4.2 and 
A4.5 of Hecht and Pang [33| adapted to our notations for [/] = [Nc] and [/] = [Nc — 1,1] 
respectively. We recall that 



C'[A^e-l,l](gU(4)) 



iVc(3iVc + 4) 
8 



(C2) 
(C3) 



From Eqs. ^-i^ and (ICTll . one can derive the following matrix elements: 
([iV, - 1, l]p; II',; SS',\Gl^\[N, - 1, l]p; II3; SS3) = 



,5+/+! 



(25+l)(2/ + i; 



'C[/1(SU(4)) 





1 




fl 1 






qi 

^3 


S'3) 


I3 
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X E {-if'-''''my[f'vmc - 1, mm'WifW'mc - 1, ib) 



xV(25^ + l) (2/^+1) 



( [/] 


[21'] 


1/1] 1 


S 1 s 


► < 




[icSc 


11 




S', 1/2 5, _ 







, (C4) 



{[N, - 1, ^^3; -^-S^iniiVc - 1, l]p; Ih- SSs) 



-l)''^'^'^\/(2'5 + l)(2/+l) 



S 1 
S3 SI 



s 



Is 



X 



E (-1)^^+'^ {KmpV'vmc - 1, 



I 1 

5 1 5 
1/2 5c 1/2 



> < 



/ 1 / 

1/2 Ic 1/2 



(C5) 



([TVe - 1, lb; //g; -^-^sl^ • G'e|[7Vc-i,i b; ^^3; = 



X 



E (-i)^=+'=i^([/1p'[/1/|[A^c - 1, ib)^([/1?'[/1?"l[A^c - 1, lb) 



xV(25^ + l)(2/^+l) 



( [/] 


[211 


1/1] 1 


'1/2 




J 1/2 


> 




11 






1/2 1 

> 




1/2 1 

> 



, (C6) 



where [/] = [Nc — 2, 1] for p = 1 and [/] = [N^ — 1] for p = 2. Sc, Ic, S'^ and are determined 
by p' ,p", q' and q" . 

There are also needed 

([AT, - 1, l]p; 11'^- SS',\s ■ Sc\[Nc - 1, lb; Ih; SSs) = Ss'^SsSr^i,i-lf^'/'^ 



X 



E (-l)'^V5e(5e + l){2Sc + 1) (i^([/V[/VI[^^c - 1, lb)) 



1/2 Sc S 
Sc 1/2 1 



,(C7) 



and 



([TV, - 1, l]p; 77^ . Tc\[Nc - 1, lb; //3; SSs) = fe5357'73(-l)'^'^'l 



X E (-l)'V^c(/c + l)(2/e + 1) (i^([/V[/V|[A^c - 1, lb))' 



P ,P 



1/2 Ic I 
Ic 1/2 1 



.(C8) 
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